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ABSTRACT 

The Local Group compact elliptical galaxy M32 hosts one of the nearest candidate 
super-massive black holes (SMBHs), which has a previously suggested X-ray counter¬ 
part. Based on sensitive observations taken with the Karl G. Jansky Very Large Array 
(VLA), we detect for the first time a compact radio source coincident with the nucleus of 
M32, which exhibits an integrated flux density of ~47.3±6.1 fiJy at 6.6 GHz. We discuss 
several possibilities for the nature of this source, favoring an origin of the long-sought 
radio emission from the central SMBH, for which we also revisit the X-ray properties 
based on recently acquired Chandra and XMM-Newton data. Our VLA observations 
also discover radio emission from three previously known optical planetary nebulae in 
the inner region of M32. 

Subject headings: galaxies: nuclei —galaxies: individual (M32)—radio continuum: 
galaxies 


1. Introduction 


It is recognized that super-massive black holes (SMBHs), commonly residing in the nuclei of 
present-day galaxies with a substantial stellar bulge, have grown the bulk of their mass during 
quasar phases, but otherwise spent the majority of their life accreting at rates well below the Ed¬ 
dington limit (e.g., Soltan 1982 : Yu &: Tremaine 2002 : Marconi et al. 2004). As such, most SMBHs 
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in the local Universe manifest themselves as low-luminosity active galactic nuclei (LLAGNs; see re¬ 
view by Ho 2008), the imprint of which is often difficult to discern. Nevertheless, dedicated surveys 
of nearby LLAGNs carried out in the past two decades have yielded high detection rates in the 


radio and X-ray bands le.g.. lNagar et al.ll2000l : IHo &: Ulvestadll200ll : iGallo et al.ll2010l : iMiller et al 


2 OI 2 I L providing important clues to their nature. The current consensus is th at LLAGNs are pow - 
ered by a radiatively inefficient, advection-dominated accretion flow lADAF: iNaravan fc Yilll994l ) 
that operate at very sub-Eddington accre tion rates, which is pro bably coupled with outflows in the 
form of jets and/or winds (see review bv lYuan &: Naravanll2014l L 


In the Local Group, strong dynamical evidence of a SMBH have been found for all of the three 
galaxies with a prominent bulge: Milky Way, M31 and M32, am ong which our own Galaxy makes 
the most compelling case (jGenzel. Eisenhauer fc Gillessenll201(ll ). Interestingly, these three nearest 
SMBHs also show the lowest Eddington ratioil] among known LLAGNs, thus supplying us with 
unique opportunities to explore the poorly understood physics of SMBHs at the most quiescent 
state. 

In the compact elliptical galaxy M32, a centr al SMBH of ^2.5 x 10^ has been inferred from 


mode ling of its circumnuclear stellar kinematics (|Verolme et al.l l2002l : Ivan den Bosch fc de Zeeuw 


2011)1 ). Additional evidence for an accreting SMBH comes from the Chandra d etection of an X-ray 


source coincident with the nucleus of M32 (|Ho. Terashima fc Ulvestadl 120031 . hereafter HTU03), 
which exhibits a power-law spectrum with a remarkably low 2-10 keV luminosity of rsj 1036 gj.g g 
corresponding to only ~10~®'® of its Eddington luminosity. On the other hand, the SMBH leaves 


stellar medium (ISM) in M32 (e.g., ] 

Jo. FilioDenko & Sareent 

1997; 

Saee. Welch & Mitchell 

1998; 

Welch & Saee 

2001 

Revnivtsev et al 

2007 

; but see 

SethIbolO). 

Even at the seemingly more promis- 


ing radio bands, where synchrotron radiation from a putative jet is expected, a firm detection in the 
literature has been absent. Gompiling Very Large Array observations taken before 2000, HTU03 
derived 5 a upper limits of 875, 30, 525 and 550 //Jy at 15, 8.4, 5 and 1.4 GHz, respectively. 

In this work, we report the detection of a compact radio source at the nucleus of M32 based 
on high-resolution, sensitive observations recently carried out with the Karl G. Jansky Very Large 
Array (VLA). We also revisit the properties of the X-ray nucleus identified by HTU03, using recently 
acquired Chandra data XMM-Newton d ata. We adopt a d istance of 780 kpc for M32, effectively 
placing it at the same distance as M31 (jYoung et al.ll2008lL Throughout this work we quote la 
uncertainties unless otherwise stated. All coordinates are given for epoch J2000. 


'^Sgr A* and M31* iiave an Eddington ratio of Lboi/LEdd ~ 10 (Baganoff et al. 2003; Li et al. 2011), where 


Z/boi is the bolo metric luminosity and I/Bdd = 1-3 x 10^° [Mbh/ALq] erg s 


the Eddington luminosity, for masses 


-i38r)i.r 1 „-i 

_ _ Bdd = r.o A rC /../.-I 

of X 10® Mq ( Ghez et al. 20081 ') and ^^ 1.4 x 10® M© (Bender et al. 2005), respectively. The present work infers 
Lboi/LBdd ~ lO"’^ ® for the SMBH in M32. 
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Observations 


2.1. Radio data 


Three VLA observations of M32 were carried out in the B-configuration on July 24, 28 and 
29, 2012 (Program ID: 12A-243; PI: L. Sjouwerman). Taken in the C-band, each observation 
has a total bandwidth of 2 GHz centered at 6.6 GHz a nd an integration time of 6 hrs. The flux 
calibrator was 3C48 with a 3% systematic uncertainty (jPerlev &: Butleiil2013l l. The observations 
were phase-referenced to the calibrator J0038+4137 every 8 minutes with a switching angle of 
1?1. The positional accuracy of our phase calibrator was <C/.T02. For each observation, the data 
was flagged, calibrated and imaged following standard procedures with the Gommon Astronomy 
Software Applicationqj (GASA, version 3.4.0). We have used the SPLIT task to average the 64 
channels of each spectral window with width = 2, and employed the GLEAN task using the Multi- 
Frequency Synthesis mode, nterms = 1, and Natural weighting to achieve maximum sensitivity. 
The resultant images were made for a primary beam of ~380^^ and sampled the sysnthesized beam 
of T.'14xl'.'09 with 0'.'3 pixels. We have also used the CONGAT task to concatenate the three 
visibility data sets to produce the final, most sensitive image (Figured^), which has an average rms 
noise of 1.2 fiJy beam“^. Notably, this is a factor of ~90 (5) lower than that achieved by HTU03 
at the neighboring frequency of 5 (8.4) GHz. 


2.2. X-ray data 

The field of M32 has been observed by Chandra with its Advanced CCD Imaging Spectrom¬ 
eter (ACIS) on nine epochs between 2000-2005. For each of these observations, the raw data was 
downloaded from the public archive and reprocessed following the standard procedure, using CIAO 
v4.3 and the corresponding calibration files. The relative astrometry among the individual ob¬ 
servations was calibrated by matching centroids of point-like sources commonly detected in the 
overlapping field of view. In addition, we have carried out deep XMM-Newton observations of M32 
(PI: Q.D. Wang) on four epochs in 2011, with a total exposure of ~400 ks. For each observation, 
the data obtained from the European Photon Imaging Camera (EPIC) was reduced following the 
standard procedure, using SAS vll.0.1 and the corresponding calibration files. Table [D gives a log 
of the X-ray data. 


^http://casa.nr ao.edu 
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3. Analysis and results 
3.1. Detection of radio sources 


Figured^ displays the inner 2b"x2b" region of M32 in the concatenated C-band image, overlaid 
by optical intensity contours (cyan) that trace the starlight. A compact source, which we designate 
Rl, is clearly present at the optical center of M32. R1 shows a peak flux density ~10.7±1.2 
^Jy beam“^, but it appears more extended than a point source. Thus we performed a single 
elliptical Gaussian fit to Rl using the CASA task IMFIT. The fit suggested that the source is 
formally larger than the synthesized beam, with a deconvolved FWHM « [(2"8 ± (y.'2) x (T.'8 ± 
(y.'2)] (~10.6 pc X 6.8 pc) and a positional angle of 130?1 ±1?9. The peak of Rl was found at 


[RA, DEC] = [00^42™'4R838 ± 0!007, +40°51'54"89 ± (y.'07], which is coincident with the 2MASS 
positio n of the nucleus, [00^42™4F825, +40°51'54'.'61], to within the uncertainty of the latter 
Ijarrett et al.ll2003l b An integrated flux density of 47.3±6.1 /iJy was measured for Rl in the 
concatenated image. The uncertainty is the quadratic sum of the error reported by IMFIT and the 
3% uncertainty in the absolute flux scale. We also measured the integrated flux density of Rl in 
the individual observations, but found no statistically significant variation. We further divided the 
full bandwidth into two 1-GHz sub-bands centering at 6 and 7 GHz to measure the spectral index, 
a, defined as cc i/“. The result was only poorly-constrained, with a ~ —2.0 ± 1.7. 


Three faint, off-nucleus sources are also evident in Figure [T^, which we designate as R2, R3 
and R4. They all lie within a projected distance of 10'' from the nucleus and have a significance 
of ~3(T in the concatenated image. Table [2] gives their positions and flux densities. Interestingly, 
all three sources are spatially coincident with pl anetary nebula (P Ne) identified in the SAURON 
integral field spectroscopic observations of M32 (ISarzi et al.ll201ll : their PNe ^5, #3 and ^7 are 
coincident with R2, R3 and R4, respectively). This is demonstrated by contrasting the white and 
cyan intensity contours in Figured^. The latter are derived from an archival Hubble Space Telescope 
(HST) WFG/F502N image (Program ID: 11714). At the position of each off-nuclear radio source, 
a compact F502N source, presumably dominated by the [O IIIJA5007 line emission, is clearly seen. 
Therefore, we suggest that the three off-nucleus sources are radio counterparts of known PNe in 
M32, which most likely arise from free-free emission (Kwok 2000). To our knowledge, this signifies 
the first detection of radio PNe at the far-side of the Local Group. The 6.6 GHz monochro¬ 
matic luminosities of the three sources, ~2 x lO^^ erg s'S are comparable to NGG7027, one of 
the most luminous Galactic PNe with a 4.9 GHz monochromatic luminosity of ~3 x 10^^ erg s~^ 
( Ziilstra. van Hoof fc Perlev 2008 ). 


= j/Li, 
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3.2. Revisiting the X-ray nucleus 


Figure [T]3 shows a 0.3-8 keV image for the same region as in Figured^, obtained by combining 
the nine Chandra/ACIS observations. Three discrete sources, which have been previously desig¬ 
nated as XI, X2 and X3 by HTU03, are clearly present, and they remain the only sources detected 
within this region, despite the ~3 times deeper total exposure achieved here. Among them, XI is 
coincident with the nucleus and has been suggested to be the X-ray counterpart of the SMBH by 
HTU03. The slight offset (~0'.'34) between the centriod of XI and the nucleus in the WFC/F502N 
image is consistent with the astrometry accurac50. From each ACIS observation, we extracted a 
spectrum for XI from a 2'hradius circle, and a background spectrum from an annulus with inner- 
to-outer radii of 2"-4!'. In the following we restrict our analysis to energies above 0.5 keV, so that 
contamination from X2, the nearby very soft source, is negligible. We fitted the individual spectra 
using an absorbed power-law model {wabs*powerlaw in XSPEC), fixing the equivalent hydrogen col¬ 
umn density at the Galactic foregro und value, 7 x 10^^ cm ~^. This is justified by the fact that M32 
exhibits no detectable neutral ISM ( Welch &: Sage 2001 '). For ObsIDs 2017, 2494 and 5690, which 


have good signal-to-noise ratios, we were able to constrain the photon-index, F, to be 2.50j’jQj7, 
2.72^Q^g and 2.36 ^q5o, respectively. For the remaining ObsIDs, we fixed the photon-index at 2.4 
and constrained the normalization only. The spectral models were then used to predict the 0.5-10 
keV intrinsic luminosity of XI at each epoch. 

The moderate angular resolution afforded by the XMM-Newton/FiPlC observations is not 
optimal for isolating the emission from XI, given the presence of the brighter source X 3 located 


9" away (Figure [TJi). Therefore, we employed a two-dimensional image fitting method (|Li et al 


20 111 ) to simultaneously determine the observed fluxes of XI and X3, which not only maximizes the 
counting statistics, but also accounts for the propagation of the dominating error arising from the 
mutual PSF scattering. Briefly, we fitted the 0.5-8 keV image of the central ~60^^ x 60^^ field with a 
three-component model convolved with the PSF. The model consists of a constant local background 
and two point sources (for XI and X3), each represented by a delta function. We mimicked the PS F 
using the image of a bright transient source recently found in the bulge of M31 (jKaur et al.ll2012l l. 
averaged over four AMM-Aewfon observations (ObsIDs 0600660301, 0600660401, 0600660501 and 
0600660601). We jointly fitted the EPIC MOSl and MOS2 images taking into account their distinct 
PSFs. We did not include data from the PN due to its poorer angular resolution. In the fit, we fixed 
the relative offset between XI and X3, as determined fr om the combined Chandra/ACIS image. 
The best-fit was obtained by minimizing the C-statistic (ICash et al.lll979l ). The resultant count 
rates and uncertainties of XI from the individual observations were then converted into an intrinsic 
luminosity by adopting the above absorbed power-law spectral model (P = 2.4). We have added 
to the error budge a 10% uncertainty in quadrature, to account for the uncertainty in the adopted 
photon-index. A consistency check indicates that the average flux ratio between XI and X3 is in 


^The typical pointing accuracy of both HST and Chandra is < 1". However, owing to the lack of common sources, 
it is virtually impossible to match the optical/X-ray astrometry to better than the individual pointing accuracy. 
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good agreement with that measured from the Chandra observations. 

The 0.5-10 keV intrinsic luminosity and its uncertainty thus determined for each observation 
are given in the last column of Tabled) Figure [2] shows the long-term light curve of XI over ^11 
yrs. The average luminosity seems to have increased by a factor of ~2 since 2002, although the 
sparse temporal sampling could have missed intervals of substantially weaker or stronger radiation. 
In addition, the luminosity was seen to increase by ~(190 ± 100)% between two XMM-Newton 
observations separated by two days (Table [T|) • We also searched for but found no significant intra¬ 
observation flux variation in either the Chandra or XMM-Newton exposures. 


4. Discussion and summary 


The VLA observations with unprecedented sensitivity have helped us establish the presence 
a compact radio source in the nucleus of M32. The source, Rl, has a monochromatic luminosity 
of ~2.3 X 10^^ erg s“^ at 6.6 GHz. This remarkably weak emission and the lack of significant 
variability hinder a straightforward identification of its nature. We discuss several possibilities 
below and favor an association with the putative SMBH. 


First, Rl can in principle arise from by synchrotron or free-free emission powered by, for 
instance, nuclear star-forming activities, in which case the source might span a physical size of 
~10 pc, as its slightly extended morphology suggests (Section l3.1ji . However, there is virtually no 
evidenc e for recent or on-g oing star formation at the center of M32 - in fact, none across the entire 
galaxy (iBrown et al.l 120011 ) . disfavoring the presence of any recent core-collasped supernovae. The 


expected rate of Type la supernovae at the nucleus is also very low (< 10“^ yr“^; HTU03). Thus, 
Rl is unlikely to be associated with synchrotron emission from a supernova remnant. Neither is 
there evidence for a circumnuclear ionized gas in M32. Ho et al. (2003) reported a 3a upper limit in 
the Ha luminosity of < 2 x 10^® erg s“^. Assuming a typical electron temperature Tg ~ 10“^ K, we 
can use the following relation to estimat e the strength of free-free emission from any circumnuclear 
ionized gas (jCaolan &: Deharveneill986l h 


F(Ha) 


10 erg cm ^ s ^ 


0.8 


(lO^Kyl 


\ - 0.52 


GHz 


xO.l 

/ ymJy 


( 1 ) 


The above upper limit in Ha luminosity corresponds to < 28.4 /rJy. This is significantly 

lower than the flux density of Rl, suggesting that the bulk of the observed radio emission does not 
come from free-free emission. 

Alternatively, we might have caught a transient radio source sim i lar to those found in the ce nter 
of our Galaxy (e.g., Davies et al. 1976 : Zhao. Roberts fc Goss 1992; Hvman et al. 2002 . 2009). In 
particular, the Galactic Genter Transient (GGT) detected at a projected distance of ~1.4 pc from 
Sgr A* during 1990-1991 would have a 6.6 GHz flux density of ~10-20 pJy if located at the distance 
of M32 (Zhao et al. 1992). The nature of some of these transients, including the GGT, remains 
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uncertain. One favored suggestion was that they are radio outbursts of black hole binaries (BHBs), 
although in most cases the expected X-ray emission was actually not detected. R1 could be such 
a radio transient, provided that it was not physically related to the stable source XI. 

Can, then, both R1 and XI be the persistent counterpart of a BHB located in the center of 
M32? Indeed, the identification of XI as the X-ray counterpart of the SMBH essentially rested on 
positional coincidence (HTU03). The very modest luminosity of XI, ~10^® erg s“^, lies within the 
range of X-ray binaries (XRBs) in the low/hard state. The Chandra and XMM-Newton observations 



log Lr = ( 4.80 ± 0 . 24 ) -k ( 0.78 ± 0 . 27 ) log Mrh + ( 0.67 ± 0 . 12 ) log Lx ( 2 ) 

Following this relation, a BHB with Mrh IOMq and L 2 _ioke 1036 

erg s is expected to 

have Lsghz of ~ 5 x 10 ^® erg s“^. This is ~400 times lower than our measurement for Rl, when 
extrapolated from the neighboring frequency of 6.6 GHz assuming a canonical spectral index of 
a = - 0 . 7 . 


Lastly, Rl may be the long-sought radio counterpart of the central SMBH, which most likely 
arises from radio synchrotron radiation. If this were the case, Rl w ould be one of the weakest 
known “AGNs” in radio, along with Sgr A* (L^qh,, 3 x 10^^ erg s“^; Zhao. Bower &: G^boOlh 


and M31* (Lsghz ~ 2 x 10^^ erg s~^: lGarcia et al.ll2010l L On the other hand, the measured value of 
r (~2.4) for XI would be compatible with the empirical anti-correlation between the photon-index 
and the Eddington ratio in LLAGNs ( Gu fc Gaa 2009 1. for Lboi/LEdd ~ 10“^'^ of the SMBH in 
M32. This supports XI being the genuine X-ray counterpart of the SMBH. 


Now, assuming that both Rl and XI are the manifestations of the central SMBH, it would 
be interesting to confront the radio and X-ray measurements with the FP. The “classical” FP 
(Equation 2) predicts Lr ~ 8.1 x 10^^ erg s“^ (for Mrh ~ 2.5 x 10® Mq and Lx ~ 1.0 x 10®® 
erg s“^). Alternatively, the modified EP proposed by Yuan & Gui (2005) and Yuan et al. (2009; 
Equation 5 therein). 


logLR = 1.29(±0.03)logLx + 0.11(±0.04)logMBH - 14.1, (3) 

which should operate at the low-luminosity regime (Lx < 10“®LEdd) where the X-ray emission 
might be dominated by a synchrotron-cooled jet, predicts Lr j'NJ 1 .1 X 10®® erg s Taking into 


^Scaled from the average 0.5-10 keV luminosity assuming the fiducial power-law spectrum iSection 13.211 . 































account the scatter (~1.0 dex; iGiiltekin et alJ l2009|) in the FPs, the predicted values are not 
inconsistent with the observed luminosity (2.3 x 10^^ erg s~^). This can be further considered a 
supporting evidence that both R1 and XI are indeed physically related to the central SMBH. In 
this case, the slightly extended morphology of R1 might trace pc-scale jets/lobes or a core-plus-jets 
system. 


In summary, our sensitive VLA observations have revealed a weak radio source coincident with 
the nucleus of M32. While at this stage a stellar origin cannot be completely ruled out for this 
nuclear source, we favor the interpretation that it is the long-sought radio counterpart of the central 
SMBH, and suggest to designate it M32* for future reference, following the convention of Sgr A*, the 
radio counterpart of the Galactic Genter Black Hole. Future VLA observations in multi-wavelength 
and with higher angular resolutions will be useful to unambiguously determine the nature of this 
source. Such observations will also benefit an in-depth investigation of the serendipitously detected 
radio PNe in the vicinity of M32*. 
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Fig. 1.— (a) The 18-hr integrated VLA C-band image of the central 25" x 25" (~95 pc x 95 
pc) region of M32, overlaid by intensity contours (positive in white solid and negative in blue 
dashed) at levels of —3, —2, 4, 5, 7 and 8 times the average rms (1.2 /rJy beam“^). Also shown are 
HST/WFC3/F502N intensity contours (cyan) highlighting the starlight distribution of M32. Four 
compact sources are labeled, among which R1 is coincident with the nucleus of M32 (named M32*), 
while R2, R3 and R4 are coincident with optical PNe. The black ellipse on the lower left-hand 
corner represents the synthesized beam size, (b) A Chandra 0.3-8 keV image showing the same 
region, overlaid by the F502N intensity contours. The three X-ray sources are labeled, following 
the convention of HTU03. 
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Fig. 2.— A long-term 0.5-10 keV light curve of XI based on nine Chandra (triangles) and four 
XMM-Newton (squares) observations. See text for details. 


Table 1. X-ray observations of M32 


ObsID 

Date 

Exposure 


(1) 

(2) 

(3) 

(6) 

313 (C/ACIS-S) 

2000-09-21 

6.0 

^■0-oi 

314 (C/ACIS-S) 

2000-10-21 

5.1 

1 0+0.4 
-*-•^-0.4 

1580 (C/ACIS-S) 

2000-11-17 

5.1 

-‘-•0-0.4 

1584 (C/ACIS-I) 

2001-07-03 

4.9 

1 Q+U.d 

-*-•^-0.3 

2017 (C/ACIS-S) 

2001-07-24 

45.9 

1 ’Y+U.'J 

‘ -0.2 

2494 (C/ACIS-S) 

2001-07-28 

16.0 

1 y+O.S 
-*-•'-0.3 

1576 (C/ACIS-I) 

2001-10-05 

4.9 

1 O+0.4 
-*-•^-0.4 

2894 (C/ACIS-I) 

2002-10-14 

4.7 

^•^-0.8 

5690 (C/ACIS-S) 

2005-05-27 

113.0 

0 4+u-^ 
^•^-0.3 

0672130101 (X) 

2011-06-27 

99.4 

1 Q+0-3 
-*-•^-0.3 

0672130601 (X) 

2011-07-05 

119.9 

1 <7+0.3 

-*-•^-0.3 

0672130701 (X) 

2011-07-07 

101.7 

Q O+0.4 
'^•^-0.4 

0672130501 (X) 

2011-07-13 

55.8 

ry O+0.4 
^•®-0.4 


Note. — (1) Observation ID. Chandra observations 
are denoted by “C/ACIS-S” or “C/ACIS-I” and XMM- 
Newton observations by “X”; (2) Date of observation; (3) 
Effective exposure, in ks; (6) 0.5-10 keV unabsorbed lu¬ 
minosity of the nuclear source, in units of 10^® erg s“^. 
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Table 2. Properties of detected radio sources 


Source 

(1) 

RA 

(2) 

DEC 

(3) 

-^P 

(4) 

56.6 

(5) 

Le.e 

(6) 

l'[OIII] 

(7) 

R1 

00'‘42’"4H838 

+40°51'54'.'98 

10.7±1.2 

47.3±6.1 

22.7 

< 4.3 

R2 

00'‘42"*4R891 

+40°51'58'.'41 

6 .6±1.2 

3.5±1.2 

1.7 

11.1 

R3 

OO'* 42"* 42; 563 

+40°51'56'.'52 

6 .6±1.2 

3.4±1.2 

1.6 

7.0 

R4 

OO'* 42"* 42; 659 

+40°51'56'.'41 

6 .8±1.2 

3.5±1.2 

1.7 

5.6 


Note. — (1) Source name; (2)-(3) Source position in Celestial coordinates; (4) Peak 
flux density, in units of pjy beam“^; (5) Total integrated flux density and error at 6.6 
GHz, in units of jJ,3y, (6) Monochromatic luminosity at 6.6 GHz, in units of 10®^ erg s“^; 
(7) [O HI]A5007 emission line luminosity (or 3a upper limit), in units of 10®® erg s“^, 
from Sarzi et al. (2011). 






